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The influence of next-nearest neighbor electron hopping, t', on the polaron 
and bipolaron formation in a square Hubbard-Holstein model is investigated 
within a variational approach. The results for electron-phonon and electron- 
electron correlation functions show that a negative value of t' induces a strong 
anisotropy in the lattice distortions favoring the formation of nearest neigh- 
bor intersite bipolaron. The role of t' , electron-phonon and electron-electron 
interactions is briefly discussed in view of the formation of charged striped 
domains. 
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I. INTRODUCTION 



In recent years the experimental evidence in favor of polaronic carriers in doped cuprates 
and in manganese oxide perovskites has grown. In manganites a large amount of experi- 
mental results, ranging from EXAFS0 measurements of lattice distortions to giant isotope 
shift of the Curie temperaturei and to frequency shifts of the internal phonon modesjl 
have pointed out the relevance of the Jahn- Teller polaron formation^ beside the double and 
super-exchange magnetic effects.il Also in cuprates there is strong experimental evidence 
supporting a relevant role of the interaction between charge carriers and lattice distortions 
in addition to the strong electron correlations. Optical experiments in the mid-infrared fre- 
quency region,! atomic pair distribution function analysis of the neutron powder diffraction 
datJ and of the EXAFS signals! due to the Cu — O bond distances have shown polaronic 
effects in doped cuprates pointing out a strong response of the local structure to the charge 
state. 

This large amount of experimental data has renewed the interest in studying problems 
of electrons interacting with the lattice degrees of freedom. In literature several models 
have been introduced to treat the electron-electron (el-el) and electron-phonon (e-ph) in- 
teractions in these compounds.i'0 In this paper we will restrict our attention to one of 
the most simple and frequently considered models for the polaron and bipolaron formation: 
the Holstein-Hubbard model. The problem of a single tight binding electron coupled to an 
optical local phonon mode has been analyzed in several ways: Monte Carlo simulations,0 
numerical exact diagonalization of small clusters,!!! dynamical mean field theoryjll density 
matrix renormalization groupEl and variational approaches. EUll As a common result the 
ground state energy and the effective mass in the Holstein model are continuous functions 
of the e-ph coupling constant and this one-body system does not have phase transitional 
Depending on the adiabatic parameter the ground state properties can change more or less 
significantly but without breaking the translational symmetry. Recently the influence of 
the Hubbard repulsion on bipolaron formation has been investigated^'^ by variational and 
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exact diagonalization methods in one and two dimensions. In the adiabatic regime the pos- 
sibility of formation of intersite bipolarons has been suggested and it has been shown that 
their mass is significantly reduced with respect to the on-site bipolaron. 

On the other hand accurate investigations on high temperature superconducting mate- 
rials have shown that some properties of cuprates compounds, as the shape of the Fermi 
surface or the band structure, can be explained introducing a next-nearest neighbor (nnn) 
electron hopping term.0 This term is essential for reproducing the experimentally observed 
behavior of the electron band near the M points of the Brillouin zone and it allows to handle 
the differences between electron and hole doped materials. The inclusion of nnn electron 
hopping is expected to affect significantly the behavior of the system. For instance, it has 
been shown that t' may deeply modify the properties of the t — j and Hubbard models 
and that the renormalization of the bare parameters, t and t', can be very strong also for 
moderate values of the electron correlations.^ 

In this paper we investigate, within a variational approach, the influence of nnn transfer 
integral on the polaron and bipolaron formation in the two dimensional Hubbard-Holstein 
model. In particular for the single electron we use a recently proposed variational approachlil 
based on a linear superposition of Bloch states that describe large and small polaron wave 
functions. This approach provides a very good description of the polaron features in any 
regime of the parameters of the Holstein model and it does not involve a truncation of the 
boson Hilbert space as required by all the numerical techniques. The computational effort 
is very little involving few variational parameters. A variational approach is used also to 
study the ground state bipolaron features. In this case we use a single wave function able to 
interpolate between large and small bipolaron regimes. In order to check the validity of the 
proposed wave function we compare our results for the one-dimensional Hubbard-Holstein 
model with those recently published by Bonca et al.,@ which are the most accurate bipolaron 
available calculations. In the most interesting regime, characterized by electron and phonon 
energy scales not well separated, the difference between the bipolaron ground state energy 
estimations of the two methods is about 0.1%: this makes us confident of the accuracy of the 
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proposed approach. Next we apply the variational method to the two dimensional extended 
Holstein-Hubbard model. 

In this paper we show that a negative value of t' induces a relevant anisotropy in the 
lattice displacements associated to the polaron formation. A decreasing of t', in the range 
of parameters of physical interest, strengthens the e-ph correlations along the (1,0) and 
(0, 1) directions and reduces the lattice displacements along the (1,1) and (1,-1) directions. 
Moreover a negative value of t' favors the mobile intersite bipolaron formation along the (1, 0) 
and (0, 1) directions and, when this does not happen, the el-el effective potential shows a 
strong dependence on the spatial directions. Although these results concern the case of only 
two electrons interacting with local phonon modes, we believe that the formation of intersite 
bipolaron in the (1,0) direction suggests the possibility that, at physical relevant doping in 
the manganites and cuprates, the e-ph interaction together with anisotropy can favor the 
formation of charged stripes. Of course, further investigation is needed to support this idea 
in specific materials. 

The paper is organized in the following way. In section 2 the model is introduced. In 
section 3 the variational approach for a single polaron is discussed. In section 4 the bipolaron 
properties within the one-dimensional Holstein-Hubbard model are analyzed by means of 
a variational method and are successful compared with the best available results recently 
published by Bonca et al..@ Successively the variational approach is extended to study 
the bipolaron features within the two-dimensional Hubbard-Holstein model. In section 5 
the numerical results are presented and the role of t', e-ph and el-el interactions is briefly 
discussed in view of the formation of charged striped domains. 

II. THE MODEL 

The two-dimensional extended Hubbard-Holstein model is described by the Hamiltonian: 

H = -t^2 C l+S,a C i,° ~ t' C i+5',a C i,<r + ^0 ^q 

i,5,tj i : 8 f ,tj q 
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+7= E 4,<V + 4*^) + f/ E ^4 • (i) 

V-'* i t<T ,q i 

In the Eq.(jl]) cf denotes the electron creator operator at site i, whose position vector is 
indicated by Ri, at represents the creation operator for phonon with wave number q, t and 
t' are, respectively, the transfer integral between nearest and next nearest neighbor sites, S 
and 5' indicate, respectively, the next and next-nearerst neighbors, u is the frequency of 
the optical local phonon mode, U represents the Hubbard interaction for electrons on the 
same site and M q indicates the e-ph matrix element. In the Holstein model (short range 
e-ph interaction) M q assumes the form: 

M q = -?=u> . (2) 

Here N is the number of lattice sites. 

III. POLARON VARIATIONAL APPROACH 

In previous papersjl! we have developed a variational approach based on a linear super- 
position of Bloch states that represent the large and small polaron wave functions. These 
two wave functions are chosen as translationally invariant Bloch states, they are obtained 
from localized states centered on different lattice sites, just like a band state is related to 
atomic orbitals, and they provide a very accurate description of the two asymptotic regimes 
of weak and strong e-ph coupling: 



l4 0) > = -7f?E^l4 0) (4)> (3) 



i 

k \/N 

V 1 » n 

where 

\^f{Rn)) = e^ )Ck)a ^ Rn - h - c ] E4 a) (^)cU m (|0) ei ® \0) ph ) . (4) 

m 

In the Eqs.(|3],£|) the apex a indicates the large (a = I) and small (a = s) polaron 
wave function, |0) e ; and \0) p h denote the electron and boson vacuum states and 0^ (i^) 
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are variational parameters such that J2 m 10? (^m)| 2 = 1 • These two wave functions are 
characterized by different phonon distribution functions: 



f!\k) = *r (5) 

u + E b {k + q)-E b {k) 



and 



4 s) w = ^E^i4 s) (^)i 2 - ( fi ) 



Here E b (q) is the free electron band energy and the variational parameters <f& (R m ) take 
into account the broadening of the electron wave function. In this paper we restrict the sum 
in Eq.(|]) till to fifth neighbors. It is evident that the large polaron wave function takes into 
account the average effect of the correlation introduced by the electron recoil (Eq||), effect 
absent in the small polaron phonon distribution function. 

We have shownll! that these wave functions, far away from the two asymptotic regimes, 
are not orthogonal and the off-diagonal matrix elements of the Holstein Hamiltonian are not 
zero. This suggests that the lowest state of the system is made of a mixture of the large and 
small polaron solutions and so justifies the idea to use a variational method to determine the 
ground state energy by considering as trial state a linear superposition of the wave functions 
describing the two types of previously discussed polarons. At f — the comparison of the 
numerical results with the data of the Density Matrix Renormalization group0 and Global 
Local Variational methods has shown the great accuracy of the proposed approach. 



IV. BIPOLARON VARIATIONAL APPROACH 

A. One-dimensional case 

Regarding the bipolaron, the most general wave function of two electrons in a periodic 
potential interacting with the longitudinal optical phonons and with each other through the 
Coulomb force can be written as: 
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_ _ .7* [Rn-i ~\- Rnn ) — » — * 

W = M £ C ^mR ni ,Rn 2 )) . (7) 

ni,n 2 

|^) is a state that is multiplied by the factor e tk ' Rm under a lattice vector R m translation 
and k is the Bloch state wave number. The Eq.(^), i.e. the polaron wave function component 
describing the charge carrier distributed around the site R n , is the key ingredient to build 
the bipolaron wave function component \i[)^(R ni , R n2 ))- Fixed the relative distance between 
the centers of the two electrons (R ni — R n2 )i we adopt a trial wave function for the singlet 
state that is proportional to the product of two polaron wave functions centered on R ni and 
R„„ sites: 



\^ k (R ni ,R n2 ))= 1% {R ni - ^^M^-^W 



-/i.e.} 



^(f)j:(R mi ,R ni - Rn 2 )c < n 1+mi ^^24>k{Rm2,Rn 1 ~ Rn 2 ) C n 2 +m 2 ,i (i )e* ® \ Q )ph) 

mi m 2 



where 0, 7 and h are variational functions. In particular, r yt(Rn 1 — R n2 ) gives the weight 
of the bipolaron wave function component with the centers of the two charge carriers at 
distance \R ni — R n2 \- For the phonon distribution function h we assume the following form: 

hJh P - f? n _ g^odjk, Rrn -R n2 )l Vn 

u* + h x {k)(E h (t,k + $-E h (t,k)) 

that is able to interpolate between the asymptotic expressions of large {h\ = d = 1) and 
small (hi = and d — 1) polaron. Here Eb(t,q) = —2tcos(q x a) is the free electron band 
energy, a is the lattice constant and d(k, R ni — R n2 ) and hi(k) are variational functions. We 
note that we have not used a linear superposition of polaron wave functions to build the 
bipolaron wave function, but, following Toyozawa,0 a single wave function able to interpolate 
between the two asymptotic regimes, as shown by Romero et al..@ This approach allows us 
to reduce the number of variational parameters, it requires, to be implemented, a very little 
computational effort and it provides a very good description of the bipolaron ground state 
properties as will be shown in the following. 

The coefficient <fi has been chosen such that it takes into account the broadening of the 
electron wave function till to third neighbors. Therefore, for any fixed relative distance 



*\Rmi R 



between the centers of the two charge carriers, we introduce three independent variational 
parameters: 

a(\R ni -R n2 \) if \R m \ =0 

P(\R ni — Rn 2 \) if \Rm\ = 0- 
- 1 l(\Rni ~ Rn 2 \) if \Rm\ = 2a 

S(\R ni - R n2 \) if \R m \ = 3a 
otherwise 

with a 2 + 2(/3 2 + 7 2 + 5 2 ) = 1. In this paper, to simplify the numerical calculations, we have 
chosen a, j3, 7, 5 independent on \R ni — R n2 \ for \R ni — R n2 \ > 3a. Likewise the function 
d(k, R ni — R n2 ) assumes the following form: 

d {k) if \R ni - R n2 \ =0 
di(k) if \R ni — R n , 2 \ = a 
d 2 (k) if |i? ni - R n2 \ = 2a 
d^(k) otherwise 

where d (k), di(k), d 2 (k) and d 3 (k) are variational parameters. Finally we choose the fol- 
lowing form for the function 7^(-R n i — Rn 2 ) 



d(k, R ni — R 



n 2 1 



IkiRnt ~ Rn 2 



if \R 



in 



Rn 2 



if \R 



in 



Rn 2 



7b(*) 

72 (k) if \R m - R n2 \ 

j 3 (k)e~' y4( - k ^ Rn ^~ R "^ otherwise 





a 

2a 



where 7o(fc), 1i(k), 72(fc) and ^(k) are variational parameters. The minimization of the 



quantity 



has been performed by making use of a routine based on a standard 



Newton algorithm. 

In Fig.l we plot the energy difference between the two-particle ground state and twice 
the one-particle ground state as function of the Hubbard repulsion U at t = ujq and g — 1. 
The comparison with the numerical results recently published by Bonca et al.@ shows the 



accuracy of the proposed variational approach. For example, the present method's estimate 
of the bipolaron ground state energy for g — 1, t — Uq and U = in the thermodynamic limit 
is Eq = —5.4185 that is in excellent agreement with the Bonca's result: Eq = —5.4246, the 
difference being about 0.1%. In Fig. lb and Fig.lc the electron-electron correlation function 

p{R m ) = (V , jf =0 l52i n * t T n i+m.llV'ifc = o)/(V'ie = olV'jfc=o) is pl° tted 101 two different values of U at 
t = ujq and g = 1. The results are compared with the data reported by Bonca et al. El. The 
comparison shows that the used wave function provides not only an excellent estimation of 
the ground state energy but gives a high accurate description of the bipolaron ground state 
properties. In the next subsection we extend the proposed approach to the two dimensional 
case. 



B. Two-dimensional case 

Here we extend the variational approach to explore the bipolaron features for electrons 
interacting with longitudinal optical phonons and with each other in a two-dimensional 
lattice. As in the one- dimensional case, we adopt a trial wave function for the singlet state 
which is a boson coherent state multiplied by the product of linear superpositions of Wannier 
wave functions: 



m{Rn 1 ,Rn 2 ))=lk{ R ni ~ Rn 2 )e " 



e i<l-Rn 1+e iq-Rn 2 _ fec _J. 

4>k{Rmn Rni ~ Rn 2 )cl 1+mi ^ X%( R ra 2 , R ni ~ Rn 2 )c ] n2+m2 ,i (l°)eZ ® \fyph) 
mi m 2 

where 0, x, 7 and h are variational functions, with Xk(Rm, R-m ~ Rn 2 ) = 4>k(Rm, Rn 2 — Rm)- 
For the phonon distribution function we assume the following form: 

hJh R _ p \ guj d(k,R ni -RnJ/y/N 

lu + h x {k) [E h {t, k + q)- E b (t, k)) + h 2 (k) (E b (t>, k + q)- E b (t>, k)) 

that is able to interpolate between the asymptotic expressions of large (hi — hi = d — 1) and 
small (hi = hi = and d = 1) polaron. Here E b (t, q) = —It (cos(q x a) + cos(q y a)), E b (t', q) = 
—At' cos(q x a) cos(q y a) and d(k,R ni — R n2 ), hi(k) and h 2 (k) are variational functions. The 
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coefficients and x have been chosen such that they take into account the broadening of 
the electron wave function till to next-nearest neighbors, i.e. the sums in the expression of 
\i[)j:(R ni , R n2 )) are restricted till to next nearest neighbors. Furthermore and x have the 

— » — * . — 

bipolaron symmetry for relative distances \R ni — R n2 \ < 2y2a (fifth neighbors), while for 
\R ni — R n2 \ > 2v^2a we assume, in analogous way to the one-dimensional case, the following 



form for 



— * — * — * 



a if \R m \ = 

P if Rm = 5 
7 if Rm = 6' 
otherwise 



with a 2 + 4(/9 2 + 7 2 ) = 1 and a, (3, 7 independent on (R ni — R n2 ). The on site, nearest 
neighbor and next-nearest neighbor values of the functions 7g(-Rm —Rn 2 ) an d h$(k, R ni —R n2 ) 
have been determined variationally whereas for relative distances between the centers of the 

— * — * . — 

two charge carriers such that \R ni — R n2 \ > y2a (next-nearest neighbors) we have used 
asymptotic expressions with parameters fixed by the variational approach as in the above- 
introduced one-dimensional case. Here we take into account the polaron and bipolaron 

— * 

ground state (k — 0). 



V. THE RESULTS 

In Fig. 2 we plot the polaron ground state energy E, the mean phonon number N, the 
spectral weight Z and the polaron kinetic energy K, in units of the bare electron kinetic 
energy, as a function of the e-ph coupling constant, at t/uo = 2, for different values of nnn 
transfer integral t'. At t' — there is a sharp transition between the large and small polaron 
solutions for A ~ 1, where A = g 2 u /4:t is the ratio between the small polaron binding 
energy and the energy gain of an itinerant electron on a rigid lattice. In the weak coupling 
regime K is ~ 1 and iV is ~ so that the electron is slightly affected by the interaction 
with the phonons: its bare mass presents weak renormalization. In the opposite regime the 

10 



polaron band collapses: the average number of phonons increases, the kinetic energy reduces 
and asymptotically tend to the values predicted by the strong coupling perturbation theory 
(N — > g 2 , K — > e~ g2 ). The ground state spectral weight shows a sharp transition at A ~ 1 
but it is evident the presence of a wide range of values of the e-ph coupling constant where 
Z is significantly smaller than the unity but not negligible: here the ground state properties 
are those characteristic of an electron weakly affected by the e-ph interaction but a large part 
of the single-particle spectral weight lies at higher energies (intermediate polaron phase@). 
Decreasing the value of t' (£' < 0), this intermediate regime becomes more important with 
respect to the small and large polaron phases: the spectral weight is equally distributed 
between ground and all the excited states and the polaron ground state features are well 
described by a linear superposition of fy^) and |V4 )■ m the weak coupling regime, for a 
fixed value of the e-ph coupling constant, decreasing t' the mean phonon number grows, the 
ground state spectral weight and the kinetic energy decrease indicating a greater polaron 
localization. The strong coupling regime shows, instead, a non monotonous behavior as 
function of t' that can be explained in terms of the third order of the strong coupling 
perturbation theory. In fact, at this order the small polaron binding energy is the sum of 
three contributions: the first term is —g 2 u)Q, the second and the third ones are proportional 
respectively to — (t 2 + t' 2 )/g 2 uo and —t 2 t'/g 3 u 2 . 

The effect of nnn transfer integral on the lattice displacements associated to the polaron 
formation, S(R m ) = (J2i n i (jk+m + a l+ m )), is particularly interesting. A negative value of 
t' induces a strong anisotropy in the lattice distortions along the (1,1) and (1,0) directions 
as it results from Fig. 3. In the range of parameters of physical interest, for a fixed value 
of the e-ph interaction, 5(1, 0)/S(0, 0) and S(l, 1)/S(0, 0) decrease in a very different way 
with t'. In fact, the ratio S(l, 1)/5(1,0) reduces in a dramatic way indicating very strong 
correlations between electron positions and lattice displacements along the symmetry axes 
of the crystal. 

A similar trend is revealed by correlation functions of two interacting electrons. We have 
investigated the bipolaron formation in the adiabatic regime where the retardation effect 
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of the e-ph interaction may favor the rise of more extended electron bound states and in 
particular the mobile intersite bipolaron formation. 00 Fig4.a shows the phase diagram for 
the transition from unbound polarons to bipolarons at fixed t and U. We stress that the 
presence of a negative nnn electron hopping extends the region of g- values where the intersite 
mobile bipolaron formation is favored. In Fig4.b the density-density correlation function 
p(Rm) = J2i n i,i n i+m, il^P) / (4>\4>) is plotted. When the intersite bipolaron formation 
takes place, the probability to find the two interacting electrons along the (1,0) and (0, 1) 
directions exceeds the 90% and the maxima of p(R m ) are located at R m = {(±1, 0), (0, ±1)}. 
When the bipolaron does not form, both along (1, 0) and (1, 1) directions p(R m ) does not 
show any structure and the typical distance between the two particles is of the order of 
the maximum allowed separation. In this case the behavior of the effective potential is 
enlightening. At the lowest order {<p{R m ) = 8 m ,o), V eff (R m ) = U5 mfi + 2^ Eg(^|(-R m ) - 
2gh^(R m )/ VN)(1 + cosq- R m ). At t' = 0, the superposition of the polarization clouds 
produces an isotropic attractive force between the two charge carriers; decreasing t', the 
effective potential is more attractive along the (1, 0) direction and less attractive along (1, 1) 
direction due to the anisotropy in the lattice distortions generated by each polaron, as it 
results from Fig4.c. 

Finally we discuss the possible consequences of our results for the formation of striped 
charge distributions in the Hubbard-Holstein model. First of all, since t' favors the formation 
of nn intersite bipolaron, we expect that the nnn electron hopping may support, for realistic 
values of the doping, the formation of striped structures arranged along characteristic di- 
rections of the crystal. Furthermore we stress the existence of an intermediate phase where 
small and large polaron properties coexist. This property of the single polaron, as discussed 
in previous papers^, can support, in a many body problem, a first order phase transition 
with coexistence of charged striped domains characterized by different densities and lattice 
deformations. The possibility of striped structures, fluctuating and short range ordered, 
and the eventuality of coexistence of two types of charge carriers forming striped domains 
have been suggested as possible scenarios for cupratesHKH. There are also strong exper- 
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imental evidences from EXAFS,ac^ neutron scatteringEH and X-ray scatterings in favor of 
the existence of charged striped structures. In both cases the interplay of t', el-ph and el-el 
interactions may significantly affect the physics of the system. Finally we want to stress that 
the Hamiltonian ([!]) represents one of the simplest models to take into account the el-el and 
e-ph interactions and that in the case of the cuprates one needs to consider more realistic 
models of electron-phonon interaction, for example, models where the Cu — O hoppings are 
modified by the ion displacements.^ However the Holstein-Hubbard model is expected to 
provide useful informations on the physics of these systems, at least on a qualitatively level, 
and it can be an interesting starting point for studying their caractheristic features. 

In conclusion we have investigated within a variational approach the influence of next- 
nearest neighbor transfer integral on the polaron and bipolaron formation in the two di- 
mensional Hubbard-Holstein model. It has been shown that a negative value of t' induces a 
relevant anisotropy in the lattice displacements associated to the polaron formation and fa- 
vors the mobile intersite bipolaron presence along the (1, 0) and (0, 1) directions. Finally we 
have discussed the relevance of these results for the formation of striped structures arranged 
along characteristic direction of the crystal and for the coexistence of charged domains char- 
acterized by different densities and lattice deformations. 
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FIGURE CAPTIONS 

Fig.l (a) The bipolaron binding energy is plotted as function of U for a one-dimensional 
lattice of 36 sites with periodic boundary conditions. The values of the parameters are 
t = u and g — 1. The stars indicate the data kindly provided by J. Bonca; (b) and 
(c): the electron-electron correlation function is plotted at g — 1 and t = u for two 
different values of U: U = in (b) and U = 1.5 in (c). The energies are given in units 
of ujq. The stars represent the data kindly provided by J. Bonca. 

Fig. 2 The polaron ground state energy, (a), the mean phonon number, (b), the spectral 
weight, (c), and the mean value of the hopping term, in units of the bare electron kinetic 
energy, (d), are reported, for a two dimensional lattice in the thermodynamic limit 
(N — > oo), at t/ujQ = 2 for different values of t'\ t'/u = (solid line), t'/uo = —0.4 
(dashed line), t'/u = —0.6 (dotted line), t'/u = —0.8 (dashed-dotted line). The 
energies are given in units of c^o- 

— * 

Fig. 3 (a): d mfi = S(\R m \)/S(\R m \ = 0) (\R m \ indicates the modulus of the vector R m in 
units of the lattice parameter a) is reported at t/uo = 2 and g = 2 for two different 
values of t'\ (b): di ;0 = S(1,1)/S(1,0) is reported at t/u = 2 for different values 
of t'\ t'/uo = (solid line), t'/uo = —0.4 (dashed line), t'/uj = —0.6 (dotted line), 
t'/uo = —0.8 (dashed-dotted line); (c): d lj0 = 5(1, 1)/S(1, 0) is reported at t/uoQ = 2 
for different values of g: g — 2 (solid line), g = 2.5 (dashed line), g = 2.65 (dotted 
line), g = 2.75 (dashed-dotted line). The data refer to two dimensional lattice in the 
thermodynamic limit (N — > oo). 

Fig. 4 (a): the phase diagram for t/uoQ = 4 and U/ujq = 30 on a square lattice 10x10 with 
periodic boundary conditions; (b): the density-density correlation function p(\R m \) is 
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plotted at t/ujQ = 4, U/u = 30 for g = 3.5 and t'/uj = —1; (c): V e ff(\R m \), measured 
by twice the one-particle potential energy, is reported at t/u = 4, U/ujq = 30, g = 2.8 
for two different values of i! . 
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